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We construct a formula for proton-nucleus total reaction cross section as a function of the mass and 
neutron excess of the target nucleus and the proton incident energy. We deduce the dependence of the 
cross section on the mass number and the proton incident energy from a simple argument involving 
the proton optical depth within the framework of a black sphere approximation of nuclei, while 
we describe the neutron excess dependence by introducing the density derivative of the symmetry 
energy, L, on the basis of a radius formula constructed from macroscopic nuclear models. We find 
that the cross section formula can reproduce the energy dependence of the cross section measured 
for stable nuclei without introducing any adjustable energy dependent parameter. We finally discuss 
whether or not the reaction cross section is affected by an extremely low density tail of the neutron 
distribution for halo nuclei. 
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I. INTRODUCTION 



The nuclear size and mass are fundamental quanti- 
ties characterizing the bulk properties of nuclei. In fact, 
the saturation of the binding energy and density de- 
duced from systematic data for masses and charge radii 
of stable nuclei reflects the behavior of the equation of 
state of nearly symmetric nuclear matter near the sat- 
uration density no. Extension of such data to a more 
neutron-rich regime opens up an opportunity of prob- 
ing the equation of state of asymmetric nuclear matter. 
One of the poorly known parameters characterizing such 
an equation of state is the density symmetry coefficient, 
L = 3nodS/dn\ no , with the nucleon density n and the 
density dependent symmetry energy S(n). The sensitiv- 
ity of L to the structure and evolution of neutron stars 
and to nuclear structure and reactions has been discussed 
in Ref. pj. This parameter is expected to be deduced 
from future systematic data for unstable nuclei because 
macroscopic nuclear models yield a clear correlation be- 
tween the root-mean-square (rms) matter radii of very 
neutron-rich nuclei and L mainly through the L depen- 
dence of the saturation density of asymmetric nuclear 
matter (2j. 

For stable nuclei, electron and proton elastic scatter- 
ing data have been obtained systematically, from which 
the charge and matter density distributions can be fairly 
well deduced 0, 0, IE IE 0' while f° r unstable nuclei, 
not only do electron elastic scattering experiments have 
yet to be performed, but also proton scattering data 
obtained from radioactive ion beams incident on a pro- 
ton target [E E| are very limited. Instead of relying on 
such elastic scattering data, deduction of the rms mat- 
ter radius of unstable nuclei has been often performed 



through measurements of reaction/interaction cross sec- 
tion (e.g., Ref. 0). However, the relation between the 
reaction/interaction cross section and the matter radius 
is not obvious theoretically and in fact is quite model 
dependent. 

Recently we systematically analyzed the proton elastic 
scattering and reaction cross section data for stable nuclei 
at proton incident energy of ~800-1000 MeV on the basis 
of a black sphere picture of nuclei 0,^3- This picture is 
originally expected to give a decent description of the re- 
action cross section for any kind of incident particle that 
tends to be attenuated in the nuclear interiors. In fact we 
showed that for proton beams incident on stable nuclei, 
the cross section of the black sphere of radius a, which 
was determined by fitting the angle of the first elastic 
diffraction peak calculated for the proton diffraction by 
a circular black disk of radius a to the measured value, is 
consistent with the measured reaction cross section [l2j . 
We also showed that fo r sta ble nuclei of mass number 
A > 50, the quantity y/3/ba agrees with the rms mat- 
ter radius deduced in previous elaborate analyses from 
the proton elastic scattering data within error bars 
These salient features of the black sphere picture could 
help deduce the density dependence of the symmetry en- 
ergy from future systematic data on the reaction cross 
section and proton elastic scattering for unstable nuclei. 

In this work we construct a formula for the proton- 
nucleus total reaction cross section as a function of the 
target mass number A, the proton incident energy T p , 
and the target neutron excess 5 = 1 — 2Z/A. The depen- 
dence of the cross section on A and T p , which is deduced 
by combining the black sphere picture of nuclei with a 
simple argument involving the proton optical depth, is 
remarkably consistent with the empirical data for stable 
nuclei at T„ = 100-1000 MeV. Bearing in mind future 
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possible data for unstable nuclei at intermediate ener- 
gies, we then describe the neutron excess dependence by 
using a radius formula Q constructed from macroscopic 
nuclear models in terms of L. As a possible application 
of the formula, we discuss the effect of nuclear deforma- 
tion and neutron halos on the reaction cross section. We 
finally extend the formula to the case of nucleus- nucleus 
reaction cross section and compare the result with the 
existing data. 



II. FORMULA FOR PROTON-NUCLEUS 
REACTION CROSS SECTION 

We start with a formula for the reaction cross section 
for protons on stable nuclei as a function of A and T„ . 
For this purpose we note the black sphere picture |ll| 
in which the black sphere radius a is determined from 
2pasin(#jv//2) = 5.1356 • • • with the proton incident mo- 
mentum in the center-of-mass frame, p, and the first peak 
angle for the measured diffraction in proton-nucleus elas- 
tic scattering, 9m- We then recall the known properties 
from the black sphere picture for T p = 800-1000 MeV 
fl2| : (i) the black sphere radius a globally behaves as 
1.214A 1 / 3 fm. (ii) the measured or agrees with no 2 
within error bars. We also note that or w tto 2 is sat- 
isfied for the values of T p down to about 100 MeV [Hil ] 
and that or is nearly flat at least for T p = 100-1000 MeV 
fl4l . It is thus natural to set 



ctr 



Aa 
a 



(1) 



where ao is the black sphere radius at T p — 800 MeV, 
and Aa is the deviation of a from ao at given T p above 
100 MeV. In this setting, we assume that the incident 
protons are point particles, leading to vanishing cross 
section for the proton-proton case. This is reasonable 
since the measured proton-proton reaction cross section 
is relatively small at T p < 1000 MeV EHH. 

In deriving the expression for Aa, we focus on the "op- 
tical" depth for incident protons. Note that in a real sit- 
uation, the black sphere radius a corresponds to a critical 
radius inside which the protons are attenuated in a tar- 
get nucleus. More explicitly, the optical depth for the 
protons, 



dl[o pn n n (r) + o pp n p (r)}, 



(2) 



amounts to a critical value when the nearest distance 
between the proton trajectory and the nuclear target is 
a. Here I is the proton trajectory, njv(r) is the point 
iV-nucleon density, and o n N is the empirical pN total 
cross section (e.g., Refs. [jjj 0, 0). Let us now ap- 
proximate the trajectory, although it is slightly distorted 
by the Coulomb repulsion, by a straight line. We also 
approximate the nucleon distributions by a trapezoidal 
form of the length of the bottom, R, and that of the 




proton 



FIG. 1: Model for the density distribution of a target nucleus 
and the critical proton trajectory inside which the reaction 
with the target nucleus occurs. 



top, R — D, as shown in Fig. 1, in such a way that they 
follow a typical behavior of the distributions deduced 
from elastic scattering data off stable nuclei [EE], i.e., 
n n {r = 0) + n p (r = 0) = p = 0.16 fm~ 3 and D = 2.2 
fm. This D corresponds to the thickness parameter of 
about 0.53 fm in the two-parameter Fermi distribution. 
We then replace the optical depth by an effective one, 



pn j 



(3) 

is the total 



T = VpNn c L', 

where o p n = (Zj A)o pp + (1 — Z/A)o 
nucleon density at r = a, and L' — 2\JR 2 — a 2 is the 
length of the part of the critical trajectory in which the 
total nucleon density is larger than n c . The value of r 
is set to be 0.9 since this is consistent with the values 
of a and n c for 12 C, 58 Ni, 124 Sn, and 208 Pb that were 
deduced from the elastic scattering data (see Ref. [Tl"| 
and references therein). 

Let us now consider the deviation AX of X = o p n, a, 
n c , and L' from the value Xq at T p = 800 MeV. As long 
as T p > 100 MeV, Ao p n is sufficiently small that 



\o p n 

T pN0 



-n c0 AL' 



L' An c 



(4) 



Because of the assumed trapezoidal distribution, one can 
obtain 



po 
D' 



Aa. 



(5) 



Combining Eqs. J2J and JSJ, one can express the relation 
between Ao p n and Aa as 



_ Ao p n 

-2 

a pN0 



po , dL' 
D L ° da 



n cQ Aa. 



(6) 
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We thus obtain 



Po«o 
Dn, 



cO 



ao dL' 
L' da 



_1 Aa 



-i 2 



pN 



(7) 



We now proceed to express the length R. In doing so, 
we first note the normalization condition for the assumed 
trapezoidal distribution, 



A = 



R 3 - 



-DR 



D 2 R - -D 3 
4 



(8) 



One can seek an approximate solution to this equation 
by setting R = R + D/2 + 5R, with R = (3^/47rp ) 1/3 , 
and assuming that 5R is small. One thus obtains 



R ~ i?o 



D 



-Roll 



12i% 
D 2 



(9) 



From this expression, 5i? can be shown to be small even 
for light elements. Hereafter we substitute expression @ 
into the reaction cross section formula JJJ through L'. 
Note that since both ao and i?o behave as oc A 1 / 3 , the 
derivative dL' /da may well be performed by assuming 
Rq ex a. 

In Fig. 2 we compare the formula Q with the empirical 
data for the reaction cross section for targets such as C, 
Ca, Zr, and Pb. The formula JJJ in which ao is set to be 
a value determined from the measured angle of the first 
diffraction maximum in proton elastic scattering [Til IT^ 
well reproduces the T„ dependence of the measured re- 



action cross section Il4| 119 for T n 



MeV. Remarkably, it agrees almost completel 
100-180 MeV data updated by Auce et al. 



down to about 100 
with the 



The T„ 



dependence of such data is stronger for smaller A, a fea- 
ture described by the prefactor of Aa P N /& p no hi Eq. J7J. 
Below 100 MeV, the measured cross section shows reso- 
nance features and sharp decrease due to the Coulomb 
repulsion between the incident proton and target nucleus, 
both of which are not allowed for in constructing the 
present formula. The neglect of the Coulomb corrections 
is reasonable above 100 MeV, since the empirical reaction 
cross section for protons on stable nuclei is appreciably 
lower than that for neutrons only at incident energy lower 
than ~ 70 MeV [H. Below 100 MeV, there is another 
limitation that faces the formula 0: a pn is too large 
to make the formula applicable |16|. We remark that 
the data above 200 MeV, tabulated in Ref. [T^. have 
to be seen with caution. This is partly because the re- 
cent data at energies between 80 and 180 MeV are 
not smoothly connected with the data above 200 MeV 
and partly because the data between 200 MeV and 600 
MeV are systematically larger than the data at 860 MeV, 
in contrast to the tendency of the nucleon-nucleon total 
cross section 0, 0] . We also remark that the previous 
predictions based on microscopic optical potential mod- 
els ^3 EH H3 are more sensitive to T p than ours in a 
way more difficult to reproduce the T p dependence of the 
data. In fact, such models involve the exponential a P N 
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FIG. 2: (Color online) The proton-nucleus reaction cross sec- 
tion as a function of proton incident energy. The target nu- 
cleus is natural C, natural Ca, 90 Zr, and natural Pb. The 
squares denote the data from Auce et al. fl9| . while the filled 
circles and crosses denote the data tabulated in Ref. [14] . Note 
some exceptions in which the target nucleus is slightly differ- 
ent in mass number: crosses for 40 Ca and natural Zr and 
squares for 208 Pb. The lines are the results from the formula 
in which a is set to be 1.2144 1 '' 3 fm (dashed lines), 2.7 
fm (bold line for C), 2.65 fm (lower, thin line for C), 2.75 fm 
(upper, thin line for C), 4.25 fm (bold line for Ca), 4.2 fm 
(lower, thin line for Ca), 4.3 fm (upper, thin line for Ca), 5.7 
fm (bold line for Zr), 5.65 fm (lower, thin line for Zr), 5.75 
fm (upper, thin line for Zr), 7.4 fm (bold line for Pb), 7.35 
fm (lower, thin line for Pb), and 7.45 fm (upper, thin line for 
Pb). The choice of 2.65-2.75 fm for C, 4.2-4.3 fm for Ca, 
5.65-5.75 fm for Zr, and 7.35-7.45 fm for Pb is based on the 
measured diffraction peak angle in proton elastic scattering. 



dependence in the prediction of the reaction cross sec- 
tion through the phase shift function, in contrast to the 
power-law a P N dependence of the formula Q ■ The im- 
portance of the power-law ct p n dependence is empirically 
known p3| . 

From Fig. 2 we can also examine to what extent a 
simple choice of ao as 1.214A 1 / 3 fm works in the cross 
section formula. As was seen from Ref. 0], this simple 
scaling well reproduces the values of ao determined from 
the measured angle of the first diffraction maximum in 
proton elastic scattering for stable nuclei ranging from 
light to heavy elements. In fact, we see a good agreement 
between the results for the cross section formula from 
both choices of ao, although the values from the simple 
A 1 / 3 scaling tend to be lower than those from the other 
by typically 5-10 % for A > 100. 

As far as the A dependence is concerned, the present 
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formula is similar to the parametrization constructed by 
Kox et al. |24(. However, they differ significantly in the 
way of describing the energy dependence. Kox et al. 
introduced a semicmpirical energy dependent parameter 
that takes care of changing "surface transparency" as the 
projectile energy changes, while the energy dependence 
of our formula stems solely from that of a P N and hence 
our formula is free from any adjustable energy dependent 
parameter. We remark that the formula by Kox et al. is 
applicable for a wider range of the incident energy includ- 
ing the low energy region where the Coulomb repulsion 
between the projectile and target takes effect, although 
it ignores resonance features which are also important in 
such a low energy region. We also remark that there is 
such a kind of formula as includes two adjustable energy 
dependent parameters |25l I26L l27j. 

We now extend the cross section formula J7J) con- 
structed for stable nuclei to unstable nuclei by introduc- 
ing the dependence on the neutron excess 5. In doing so, 
we first note a formula for the rms radii of proton and 
neutron point distributions, R p and R ni which reads 

R p = c 1 A 1 ' 3 + c 2 + c 3 {S-6 ) 2 1 (10) 

with ci = 0.915 fm, c 2 = -0.102 fm, c 3 = 0.389 fm, and 
So = 0.880, and 

R n = c 4 A^ 3 (l + c 5 LS 2 + c 6 L 2 S A ) + c 7 + c s 5, (11) 

with c 4 = 0.880 fm, c 5 = 0.00635 MeV"\ c 6 = 
-0.000172 MeV- 2 , c 7 = 0.302 fm, and c 8 = 0.193 fm. 
This formula was constructed in Ref. Q by fitting the 
calculations obtained from macroscopic nuclear models 
for A > 50, < S < 0.3, and various values of the den- 
sity symmetry coefficient L and the incompressibility Kq 
of symmetric nuclear matter. It was noted that the cal- 
culated results for the rms radii are almost independent 
of K . We remark that the formulas (|10|l and (|llf) , when 
extrapolated into the range of 4 < A < 50, reproduces 
as well the values of the rms charge and matter radii 
deduced from elastic scattering data for stable nuclei. 

The formulas (|10fl and l|ll|) imply the significance of 
allowing for uncertainties in the density symmetry coef- 
ficient L in describing the size dependent quantities at 
finite neutron excess. In fact, L acts to increase the nu- 
clear size at nonzero 6. This is mainly because the satu- 
ration density of uniform nuclear matter at finite neutron 
excess decreases from that in the symmetric case in such 
a way that the decrease is roughly proportional to LS 2 
. This effect of L could be seen by systematically an- 
alyzing the first diffraction maxima in the proton elastic 
scattering data [2|| and thus could manifest itself more 
remarkably in the reaction cross section which behaves 
as the size squared. 

In order to describe the cross section at finite neu- 
tron excess by including the effect of L, we mul- 
tiply expression JJJ) by a scale factor f(A 7 5; L) = 
R 2 n {A, 5; L) jR 2 m {A, 5 = 0; L) with the rms matter radius 

R m = ^J(N/A)Rl + (Z/A)R 2 . This is reasonable, be- 
cause, at least for stable nuclei heavier than A ~ 50, the 



rms matter radii deduced in previous elaborate analyses 
from elastic scattering data agree very well with y/Z/ha 
RlJ. Although this agreement neither holds for A < 50 
|l2| nor is obvious at large neutron excess, it is natural 
to expect that the isospin dependence of R m is similar 
to that of a. We remark in passing that the scale factor 
/ is sufficiently close to unity for stable nuclei including 
208 Pb that it does almost no damage to the good agree- 
ment between the formula and the data. 

In addition to the isospin dependence coming through 
this scale factor, another isospin dependence is inherent 
in the formula (J7J through the part associated with a P N- 
In the absence of detailed information about the surface 
diffuscness and neutron skin at large neutron excess, we 
simply leave it as it is. This is expected to be good 
enough since a pp and a pn are not so different at T p = 100- 
1000 MeV as to affect the analysis of the reaction cross 
section data within the error bars, which are typically of 
order a few percent. 

In the black sphere picture adopted here, possible neu- 
tron halos have no influence on the construction of the 
cross section formula for unstable nuclei. This is because 
such halos correspond to a very low density region which 
is "optically" very thin. Moreover, the macroscopic nu- 
clear models leading to Eqs. ifTUI) and ((TTll do not allow 
for the possible presence of a neutron halo 2] . If the cross 
section formula is significantly smaller than the empirical 
value for halo nuclei for reasonable values of L, therefore, 
this would signify the limitation of the present framework 
and the enhancement of the reaction cross section due to 
the presence of halos. 

We now focus on the interaction cross section data for 
a halo nucleus 11 Li @| as well as the reaction cross sec- 
tion data for helium isotopes [jj. Note that the interac- 
tion cross section can be regarded as equal to the reaction 
cross section for halo nuclei in which internucleon exci- 
tations generally involve breakup of halo neutrons. Such 
light elements are expected to be suitable for examining 
the isospin dependence, partly because for stable nuclei 
ao is consistent with 1.214A 1 / 3 fm, and partly because / 
can amount to about 1.1. In fact, the formula J7J with 
a Q = 1.214A 1 / 3 fm well reproduces the data for 4 He [U 
and Li of the natural abundance |3ll |. From the formula 
@ further multiplied by the scale factor / with a typical 
L value of 100 MeV, we obtain a smaller value for n Li 
by about 10 %, but a value consistent with the empirical 
data for halo nuclei 6 He and 8 He. It might be tempt- 
ing to attribute this deviation for 11 Li to the presence of 
a neutron halo, since the formula Q) and the scale fac- 
tor / arc not influenced by an extremely low density tail 
of the neutron distribution. However, this attribution is 
not obvious because of the consistency obtained for 6 He 
and 8 He which are known to have a similar three-body 
structure to that of 11 Li |3^. 

So far we have assumed that nuclei are spherical. Nu- 
clear deformation, however, does affect the reaction cross 
section and thus might be deduced from the experimental 
data measured for heavy deformed nuclei. For simplic- 
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ity, we focus on quadrupolar deformation and describe a 
nucleus as a static ellipsoid having the length, b\, 61, b 2 , 
of the three axes. Both for prolate {b\ < 62) and oblate 
(61 > 62) deformations with random orientations, the 
cross section seen by incident protons becomes, on the 
average, larger than the value for a sphere of the same 
volume. The relative increment can be calculated as 



9 = ~ 5 a 2 



(12) 



where ai is the quadrupolar deformation factor satisfying 
&2/&1 = (1 + £*2)/(l — 0*2/2). This suggests that even for 
large deformation \a%\ ~ 0.2, g amounts to about 2 % 
only. 

In the above estimate of the effect of deformation, we 
ignore possible rotation of the nucleus. Since the rotation 
speed is far smaller than the speed of light, it is reason- 
able to consider that the incident protons, which move 
relativistically, see the target nucleus nonrotating during 
the flight. However, rotation would tend to enlarge the 
nucleus itself. Balance between the centrifugal force and 
the nuclear incompressibility would typically give rise to 
about a percent enhancement of the cross section. Wc 
can thus conclude that the combined effect of deforma- 
tion and rotation is significantly smaller than the effect 
of large neutron excess. 



III. FORMULA FOR NUCLEUS-NUCLEUS 
REACTION CROSS SECTION 

In this section we generalize the formula constructed 
in Sec. [H] to the case of nucleus-nucleus reaction cross 
section. As in Ref. 01 > we straightforwardly set it as 



ctr = 7r(a P + ap) 2 , 



(13) 



where ap and ap are the black sphere radii of the projec- 
tile (Zp, Ap) and target (Zt, At) nucleus. By following 
a line of argument for the proton-nucleus case, we obtain 



<tr = 7T (ap 



Aa P )f p /2 



(aro 



Aa T )/ T /2 



(14) 



where fp and /p are the scale factors for the projectile 
and target nucleus, and Aap and Aap are the deviations 
of ap and ap from the values of apo and aro at incident 
energy per nucleon of 800 MeV, both of which are deter- 
mined from Eq. © in which a P N is replaced by 
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(15) 



for Aap and (T <-> P) for Aap. Here we have omitted the 
effects of nuclear deformation. The formula (|14fl reduces 
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FIG. 3: The 12 C- 12 C reaction cross section as a function of 
energy per nucleon. The filled circles are the empirical data 
from Refs. I24I. I30ll33|. The lines are the results from the for- 
mula 1141 in which ao is set to be 2.7 fm (bold), 2.65 fin (lower, 
thin), 2.75 fm (upper, thin), and 1.214A 1/3 fm (dashed). The 
choice of 2.65-2.75 fm is based on the measured diffraction 
peak angle in proton elastic scattering. 



to the proton- nucleus case by setting either ap or ap to 

1 /3 

be zero. Hereafter we will set apo = 1.214^" fm and 

1 /3 

aro = 1.214Aji fm unless notified otherwise. 

The formula (|14f> works well for the 12 C- 12 C case for 
energy per nucleon at least down to 200 MeV (see Fig. 
3). In fact, in this energy range the empirical values of 
the 12 C- 12 C case |24| are almost four times as large as 
those of the p- 12 C case, while the formula l|14l) for the 
12 C- 12 C case is just four times as large as that for the 
p- 12 C case. At lower energies, the formula (|14H overesti- 
mates the reaction cross section to an extent much larger 
than the one that could be explained by our neglect of the 
Coulomb repulsion. This suggests the necessity of modi- 
fying the classical optical depth argument underlying the 
formula l|14l) in analyzing the low energy nucleus- nucleus 
reactions. The agreement for energy per nucleon down 
to 200 MeV and the deviation for lower energy can also 
be seen in all the other reactions between stable nuclei 
tabulated in Table II of Ref. 0|. 

We note that for the 4 He- 12 C case, the formula (fTl|) 
significantly overestimates the reaction cross section even 
at energy per nucleon of order or higher than 800 MeV 
as can be seen by comparison with the empirical data 
[33 ■ This exceptional behavior is attributable to the fact 
that excitations associated with internucleon motion are 
highly suppressed in a particles. 

It is important to compare the constructed formula 
with the empirical data for unstable nuclei incident on a 
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FIG. 4: (Color online) Deviation of the empirical reaction 
cross section from the cross section formula 1141 with L = 
MeV (crosses) and L = 100 MeV (squares), apo and aro are 



set to be 1.2UAp /3 fm and 1.214AJ,' fm, respectively. The 
dotted lines are drawn to guide the eye. 
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stable nucleus target. We note the data for the interac- 
tion cross section at energy per nucleon of order 800 MeV 
tabulated in Ref. 01 ■ For candidates of halo nuclei such 
as 6 He, 8 He, 11 Li, 11 Be, and 19 C, we may assume that the 
interaction cross section is almost equal to the total re- 
action cross section. For a typical value of L of order 100 
MeV, we obtain a good agreement between the formula 
and the data for Be, C, and Al targets; the difference is 
of order 5 %, typically. The results from these compar- 
isons, together with the aforementioned proton-nucleus 
cases ( n Lipi, 4 He Hip, 6 He @, 8 He 0, natural Li 
|3lj . 9 Be |3lj|. natural C |3l|1. are summarized in Fig. 4. 
We find that the deviation between the empirical data 
and the formula with L = 100 MeV is within ±10 %, 
while the formula with L = MeV tends to underesti- 
mate the cross section at large neutron excess. Judging 
from the scattering pattern in the plot for L — 100 MeV, 



which is almost uniform within the band between ±10 
%, and recalling that the formula is not influenced by 
the presence of neutron halos, we see no absolute need 
for taking into account the possible enhancement of the 
reaction cross section by an extremely low density tail 
of the neutron distribution in explaining the empirical 
values of the reaction cross section. In the absence of 
this enhancement, it is expected that one can deduce the 
value of L from the comparison between the empirical 
data and the formula. We note, however, that the differ- 
ence of the formula with L = 100 MeV from the formula 
with L — MeV is of order or even smaller than that 
from the empirical data. This suggests the necessity of a 
systematic analysis in deducing the value of L itself. 



IV. CONCLUSION 

We constructed a formula for the proton-nucleus reac- 
tion cross section in a way free from any adjustable T p de- 
pendent parameter. The dependence of the cross section 
on A and T p was deduced by combining the black sphere 
picture of nuclei with a simple argument of the proton 
optical depth. For stable nuclei, this formula remarkably 
well reproduces the empirical T p dependence of the reac- 
tion cross section at T p = 100-1000 MeV. This suggests 
a great advantage over previous microscopic optical po- 
tential models in which the T p dependence becomes too 
strong through the phase shift function to reproduce the 
empirical behavior. In addition, we took into account the 
effects of large neutron excess typical of unstable nuclei 
and nuclear deformation on the reaction cross section. 
We thus found out a possible way of examining the influ- 
ence of the density dependence of the symmetry energy 
and the presence of neutron halos on the reaction cross 
section. 
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